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ABSTJUVCT 

An  analytical  method  has  been  developed 
to  determine  gear  tooth  surface  coordinates  of 
face-milled  spiral  bevel  gears.  The  method  com¬ 
bines  the  basic  gear  design  parameters  with  the 
kinematical  aspects  for  spiral  bevel  gear  manu¬ 
facturing.  A  computer  iprogram  was  developed  to 
calculate  the  surface  coordinates.  From  this  data 
a  three-dimensional  model  for  finite  element  ana¬ 
lysis  can  be  determined.  Development  of  the 
modeling  method  and  an  e.-cample  case  are  presented. 

INTRODUCTION 

Spiral  bevel  gears  are  currently  used  in  all 
helicopter  power  transmission  systems.  This  type 
of  gear  is  required  to  turn  the  corner  from  a 
horizontal  engine  to  the  vertical  rotor  shaft. 
These  gears  carry  large  loads  and  operate  at  high 
rotational  speeds.  Recent  research  has  focused  on 
understanding  many  aspects  of  spiral  bevel  gear 
operation,  including  gear  geometry  fl-121 gear 
dynamics  f 13-151 ,  lubrication  fl61 .  stress  analy¬ 
sis  and  measurement  fl7-211 .  misalignment  122.231 . 
and  coordinate  measurements  [24,251 ,  as  well  as 
other  areas. 

Research  in  gear  geometry  has  concentrated  on 
understanding  the  meshing  action  of  spiral  bevel 
gears  [8-111 .  This  meshing  action  often  results 
in  much  vibration  and  noise  due  to  an  inherent 
lack  of  conjugation.  Vibration  studies  [261  have 
shown  that  in  the  frequency  spectrum  of  an  entire 
heli-opter  transmission,  the  highest  response  can 
be  that  from  the  spiral  bevel  gear  mesh.  There¬ 
fore  if  noise  reduction  techniques  are  to  be  im¬ 
plemented  effectively,  the  meshing  action  of 
spiral  bevel  gears  must  be  understood. 

Also,  investigators  [18,191  have  found  that 
typical  design  stress  indices  for  spiral  bevel 
gears  can  be  significantly  different  from  those 
measured  experimentally.  In  addition  to  making 
the  design  process  one  of  trial  and  error  (forcing 
one  to  rely  on  past  experience),  this  inconsis¬ 
tency  makes  extrapolating  over  a  wide  range  of 
sizes  difficult,  and  an  overly  conservative  design 
can  result. 

The  objective  of  the  research  reported  herein 
was  to  develop  a  method  for  calculating  spiral 


bevel  gear-tooth  surface  coordinates  and  a  three- 
dimensional  model  for  finite  element  analysis. 
Accomplishment  of  this  task  requires  a  basic  un¬ 
derstanding  of  the  gear  manufacturing  process, 
which  IS  described  herein  by  use  of  differential 
geometry  techniques  [11 .  Both  the  manufacturing 
machine  settings  and  the  basic  gear  design  data 
were  used  in  a  numerical  analysis  procedure  chat 
yielded  the  tooth  surface  coordinates.  After  the 
tooth  surfaces  (drive  and  coast  sides)  were  de¬ 
scribed,  a  three-dimensional  model  for  the  tooth 
was  assembled.  A  computer  program  was  developed 
to  automate  the  calculation  of  the  tooth  surface 
coordinates,  and  hence,  the  coordinate  for  the 
gear-tooth  three-dimensional  finite  element  mouei. 
The  basic  development  of  the  analytical  model  is 
explained,  and  an  example  of  the  finite  element 
method  is  presented. 

DETERMINATION  OF  TOOTH  SURFACE  COORDINATES 

The  spiral  gear  machining  process  described 
in  this  paper  is  that  of  the  face-milled  type. 
Spiral  bevel  gears  manufactured  in  this  way  are 
used  extensively  in  aerospace  power  transmissions 
(i.e.,  helicopter  main/tail  rotor  transmissions) 
to  transmit  power  between  horizontal  gas  turbine 
engines  and  the  vertical  rotor  shaft.  Because 
spiral  bevel  gears  can  accommodate  various  shaft 
orientations,  they  allow  greater  freedom  for  over¬ 
all  aircraft  layout. 

In  the  following  sections  the  method  of  de¬ 
termining  gear-tooth  surface  coordinates  will  be 
described.  The  manufacturing  process  must  first 
be  understood  and  then  analytically  described. 
Equations  must  be  developed  that  relate  machine 
and  workpiece  motions  and  settings  with  the  basic 
gear  design  data.  The  simultaneous  solution  of 
these  equations  must  be  done  numerically  since  no 
closed-form  solution  exists.  A  description  of 
this  procedure  follows. 

Gear  Manufacture 

Spiral  bevel  gears  are  manufactured  on  a 
machine  like  the  one  shown  in  Fig.  1.  This 
machine  cuts  away  the  material  between  the  concave 
and  convex  tooth  surfaces  of  adjacent  teeth  simul¬ 
taneously.  The  machining  process  is  better  illus¬ 
trated  in  Fig.  2.  The  head  cutter  (holding  the 
cutting  blades  or  the  grinding  wheel)  rotates 
about  its  own  axis  at  the  proper  cutting  speed. 
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Figure  2.— Oneniation  ol  workpiece  to  generation  macPcer/. 


about  its  own  axis  at  the  proper  cutting  speed, 
independent  of  the  cradle  or  workpiece  rotation. 
The  head  cutter  is  connected  to  the  cradle  through 
an  eccentric  that  allows  adjustment  of  the  axial 
distance  between  the  cutter  center  and  cradle 
(machine)  center,  and  adjustment  of  the  angular 
position  between  the  two  axes  to  provide  the  de¬ 
sired  mean  spiral  angle.  The  cradle  and  workpiece 
are  connected  through  a  system  of  gears  and 
shafts,  which  controls  the  ratio  of  rotational 
motion  between  the  two  (ratio  of  roll).  For 
cutting,  the  ratio  is  constant,  but  for  grinding, 
it  is  a  variable. 

Computer  numerical  controlled  (CMC)  versions 
of  the  cutting  and  grinding  manufacturing  pro¬ 
cesses  are  currently  being  developed.  The  basic 
kinematics,  however,  are  still  .maintained  for  the 
generation  process;  this  is  accomplished  by  the 
CHC  machinery  duplicating  the  generating  motion 
through  point-to-point  control  of  the  machining 
surface  and  location  of  the  workpiece. 

Coordinate  Transformations 

The  surface  of  a  generated  gear  is  an 
envelope  to  the  family  of  surfaces  of  the  head 
cutter.  In  simple  terms  this  means  that  the 
points  on  the  generated  tooth  surface  are  points 
of  tangency  to  the  cutter  surface  during  manu¬ 
facture.  The  conditions  necessary  for  envelope 
existence  are  given  kinematically  by  the  equation 
of  meshing.  This  equation  can  be  stated  as  fol¬ 
lows:  the  normal  of  the  generating  surface  must 
be  perpendicular  to  the  relative  velocity  between 
the  cutter  and  the  gear-tooth  surface  at  the  point 
in  question  fll . 

The  coordinate  transformation  procedure  that 
will  now  be  described  is  required  to  locate  any 
point  from  the  head  cutter  into  a  coordinate 
system  rigidly  attached  to  the  gear  being  manu¬ 
factured.  Homogeneous  coordinates  are  used  to 


allow  rotation  and  translation  of  vectors  simply 
by  multiplying  the  matrix  transformations.  The’ 
method  used  for  the  coordinate  transfor.mation  ca 
be  found  in  Refs.  1,  S,  8  to  11,  and  27. 

To  begin  the  process  we  start  witn  the  head 
cutter  (Fig.  (8)).  In  this  report  it  is  asstuned 
that  the  cutters  are  straight-sided.  The  para¬ 
meters  u  and  9  determine  the  location  of  a 
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current  point  in  the  cutter  coordinate  system  S  . 
Angles  and  are  the  blades  that  cut  the*"" 

convex  and  concave  sides  of  the  gear  tooth  respec¬ 
tively.  Thus  r^  is  given  by  the  following 
equation: 


root  ]5  -  u  cos  jS 
u  sinsf  sin  5 
u  simf  cos  8 
1 


(1) 


Once  we  have  r^  we  then  transform  from  one  coor¬ 
dinate  system  to  the  next  for  the  coordinate  sys¬ 
tems  as  shown  in  Figs.  (4)  and  (5)  f271 . 


Using  matrix  transformations  we  can  determine 
the  coordinates  in  of  a  point  on  the  generat¬ 
ing  surface  by: 

Here  (Mjj)  describes  the  required  homogeneous 
coordinate  transformation  from  system  “j"  to 
system  "i".  Therefore  Eq.  (2)  describes  the 
location  of  a  point  in  the  gear  fixed  coordinate 
system  based  oh  machine  settings  q,  s.  r, 

5>,  and  gear  design  information  jt  and  6.  At 
this  point  the  machine  settings  and  the  gear 
design  values  are  known.  Parameters  u,  6,  and 
are  the  unknown  variables  that  are  solved  for 
numerically. 


Figure  5.— Top  and  (row  view  ol  right  hand  gear  surface  generarcr) 
(Oc  •  0  shown  here). 


Tooth  Surface  Coordinate  Solution  Procedure 

In  order  to  solve  for  the  coordinates  of  a 
spiral  bevel  gear-tooth  surface,  the  following 
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Figure  4,— Orieniaiion  of  cutler,  cradle,  and  fixed  coordinate 
systems.  Sc.  Sj.  and  Sm  respectively. 


items  must  be  used  simultaneously:  the  trans¬ 
formation  process,  the  equation  of  meshing,  and 
the  basic  gear  design  information.  The  trans¬ 
formation  process  described  previously  is  used  to 
determine  the  location  of  a  point  on  the  head  cut¬ 
ter  in  coordinate  system  S^.  Since  there  are 
three  unknown  quantities  (u.  6,  and  (5^).  three 
equations  relating  them  must  be  developed. 

Values  for  u.  9,  and  are  used  to  satisfy 
the  equation  of  meshing  given  by  Refs.  1  and  9. 

n  •  V  a  o  (3) 

where  n  is  the  normal  vector  to  the  cutter  and 
workpiece  surfa£es  at  the  specified  location  of 
interest,  and  V  is  the  relative  velocity  between 
the  cutter  and  workpiece  surfaces  at  the  specified 
location. 

Gear  design  information  is  then  used  to 
establish  an  allowable  range  of  values  of  the 
radial  (r)  and  axial  (z)  positions  that  are  known 
to  exist  on  the  gear  being  generated.  This  is 
shown  in  Fig.  6. 

First  the  equation  of  meshing  must  be  satis¬ 
fied.  This  is  given  as  f91 : 

(u  -  rcot5!cos5S)cos7sin^5  ^  q  *  p^'j 

*  ■  3in7|co3  5>sin0  r  cos7sin^sin^q  -  p^ 

t  E^|coo7oin^  +  oinTcoapcosld  -  q  ±  p^'^ 

-  I,g|oin7coo^oin^5  -  q  i  p^  "  0 
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Figure  6.— Ofientation  of  gear  le  be  gcnemicd  wild  assumed  posilions 
I  and  z. 

The  upper  and  lower  si;}ns  preceding  the 
above  terms  pertain  to  left  and  right  hand  gears 
respectively. 

The  axial  position  must  match  the  value  found 
from  transforming  the  cutter  coordinates  to 

workpiece  coordinates  S^.  fhis  is  satisfied  by 
the  following  (Fig.  6); 

z  -  E  .0  «) 

Finally  the  radial  location  from  the  work 
axis  of  rotsrion  must  be  satisfied.  This  is 
accomplishe'..  by  using  the  magnitude  of  the  loca¬ 
tion  in  question  in  the  x,  -  y^  plane  (Fig.  6): 


How  a  system  of  three  equations  (Eqs.  (3,  to 
(5))  is  solved  simultaneously  for  the  three  para¬ 
meters  u,  9,  and  tor  a  given  gear  design  with 
a  set  of  machine  tool  settings.  These  are  non¬ 
linear  algebraic  equations  that  can  be  solved 
numerically  with  commercially  available  mathe¬ 
matical  subroutines.  These  equations  are  then 
solved  simultaneously  for  each  location  of  in¬ 
terest  along  the  tooth  flank,  as  shown  in  Fig.  7. 
From  the  surface  grids,  the  active  profile 
(working  depth)  occupied  by  a  single  tooth  is 
defined. 


Figutc  7.— Calculation  points  (too  each  side)  lor  concave  and 
convex  sides  ol  tooth  surlacc. 


In  sufivmary  the  procedur  is  as  follows. 

Known  locf'ions  f  and  £  ca  the  active  profile 
are  used  along  with  the  equat„oh  of  meshing  to 
determine  unknown  parameters  b,  d,  and  <j>  .  The 
parameter  values,  machine  tool  -ettings,  and  gear 
design  values  are  used  in  the  coirdinate  trans¬ 
formation  shown  earlier  to  find  »..  e  radial  and 
axial  positions  in  the  gear  fixed  i  :?nrdinate 
system  This  procedure  results  in  the  solu¬ 

tion  of  three  simultaneous  nonlinear  algebraic 
equations  that  are  solved  numerically. 

APPLICATION  OF  SOLUTION  TECHNIO.UE 

An  application  of  the  techniques  previously 
discussed  will  now  be  presented.  The  component  to 
be  modeled  was  from  the  NASA  Lewis  Spiral  Bevel 
Gear  Test  Facility.  A  photograph  of  the  spiral 
bevel  gear  mesh  is  shown  in  Fig.  8,  and  the  design 
data  for  the  pinion  member  are  shown  in  Table  1. 


Figure  8.— Phoiograph  ol  NASA  spiral  bevel 
gear  lesi  rig  componcr.is. 

Surface  Coordinate  Calculation 

Using  Figs.  7  and  9  as  references,  we  will 
describe  the  calculation  procedure  for  surface 
coordinates.  First,  the  concave  side  of  the  tooth 
is  completely  defined  before  moving  to  the  convex 
side.  These  points  are  calculated  by  starting  at 
the  toe  end  and  at  the  lowest  point  of  active 
profile.  Nine  steps  of  equal  distance  ace  used 
from  the  beginning  of  the  active  profile  to  the 
face  angle  (addendiun)  of  the  gear  tooth,  and  then 
back  to  the  next  axial  position  (Fig,  9),  The 
procedure  is  repeated  until  the  concave  side  is 
completely  described.  Then  the  same  procedure  is 
followed  for  the  convex  side. 

Example  Model  and  Results 

From  the  one-tooth  model  described  earlier 
the  analysis  techniques  can  be  demonstrated.  The 


TABLE  I.  -  EXAMPLE  CASE  FOR  SURFACE  COORDIHATE  GEHERATIOH 


Gear  design  data 


-Number  of  teech  (pinion,  gear)  .  12.  36 

iDedendum  angle,  6,  deg  .............  .  1.0 

Addendum  angle,  deg  . 3.883 

Pitch  angle,  /i.  deg  .  18.433 

; Shaft  angle,  deg .  90.0 

Mean  spiral  angle,  deg .  35. 0 

Face  width,  mm.  (in.)  .  25.4  (1.0) 

Mean  cune  distance,  irm  (in.)  . 31.05  (3.191) 

Inside  radius  of  gear  blank,  not  (in.)  .  15.3  (0.6094) 

Top  land  thickness,  ntn  (in.) .  2.032  (0.080) 

Clearance  nsn  (in.)  . . .  6.762  (0.030) 


Generation  machine  settings 


Blade  angle.  deg 
Vector  su.m.  L^^.  .t.t..  in. 
Machine  offset.  E^^.  mm  (in.) 
Cradle  to  cutter  distance. 

s.  mm  (in.) 

Cradle  angle,  q.  deg 


Concave 

Convex 

75.222  (2.9615) 

78.1329  (3.0751) 

161.358 

24.932 

1.0363  (0.0408) 

-1.424S  (-,'.0561) 

3.9802  (0.1567) 

-ii.4856  (-V.1765* 

74.839  -2.9646) 

71.247  (2  3050. 

64.01 

53.32 

model  shown  in  Fig.  10  is  that  for  a  constant 
fillet  and  root  radius  (also  called  full  fiHet) 
model.  The  fillet  and  root  radius  on  the  convex 
side  has  been  added  along  with  the  tooth  section 
(without  the  tooth)  to  make  the  model  symmetric 
about  the  tooth  centerline.  Figure  10  shows  a 
hidden  line  plot  of  the  finite  element  meshrwith 
eight-noded  isoperimetric  three-dimensional  solid 
continuum  elements.  This  model  has  765  elements 
and  IIZO  nodes.  The  boundary  conditions  are  also 


shown  in  Fig.  10.  A  i724-HPa  (250  ksi)  constant 
pressure  load  was  applied  normal  to  the  tooth 
surface  of  nine  elements,  and  the  two  edge  sur¬ 
faces  of  the  gear  rim  had  all  degrees  of  freedom 
constrained. 

The  results  were  calculated  by  MSC/tlASTRAM 
and  were  subsequently  displayed  by  P.ATRAH.  Fig¬ 
ure  11  shows  the  major  principle  stresses  tor  the 
boundary  conditions  shown  in  Fig.  10. 
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Figure  1 1.— Major  principle  siress  for  me  boundary  conditions  spedSed  in  figure  lO. 


Sm-iMARy  OF  RESULTS 

A  method  has  been  presented  that  uses  dif¬ 
ferential  geometry  techniques  to  calculate  the 
surface  coordinates  of  face-milled  spiral  bevel 
gear  teeth.  The  coordinates  must  be  solved  for 
numerically  by  a  simultaneous  solution  of  non- 
lin‘-ar  algebraic  equations.  These  equations 
relate  the  kinematics  of  manufacture  to  the  gear 
design  parameters.  Coordinates  for  a  grid  of 
points  are  determined  for  both  the  concave  and 
convex  sides  of  the  gear  tooth.  These  coordinates 
are  Chen  combined  to  form  the  enclosed  surface  of 
one  gear  tooth.  A  computer  program,  was  developed 
to  solve  for  the  gear-tooth  surface  coordinates 
and  provide  input  to  a  three-dimensional  geometric 
modeling  program  .i.e.,  PATRAII).  This  enables  an 
analysis  by  the  finite  element  method.  An  example 
of  the  technique  was  presented. 
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